Noncollinear magnetism and spin-orbit coupling in 5d pyrochlore oxide Cd20s207 



<N 

O 

<N 

in 



U 
> 

O 

(N 
> 



Hiroshi Shinaoka,^'^ Takashi Miyake/'^ and Shoji Ishibashi^'^ 

^ Nanosystem Research Institute "RICS", National Institute of Advanced 
Industrial Science and Technology (AIST), Umezono, Tsukuba 305-8568, Japan 
^ Japan Science and Technology Agency (JST), CREST, Honcho, Kawaguchi, Saitama 332-0012, Japan 

(Dated: April 17, 2012) 

We investigated the electronic and magnetic properties of the pyrochlore oxide Cd2 0s207 using 
the density- functional theory plus on-site repulsion {U) method, and depict the ground-state phase 
diagram with respect to U . We conclude that the all- in/ all- out non-collinear magnetic order is stable 
in a wide range of U . We also show that the easy-axis anisotropy arising from the spin-orbit (SO) 
coupling plays a significant role in stabilizing the all-in/all-out magnetic order. A pseudo gap was 
observed near the transition between the antiferromagnetic metallic and insulating phases. Finally, 
we discuss possible origins of the peculiar low-temperature (T) properties observed in experiments. 

PACS numbers: 75.47.Lx,71.15.Mb,71.70.Ej 



Pyrochlore transition- metal oxides A2B2O7 are a class 
of materials that have been extensively investigated for 
several decades [1]. Recently, particular attention has 
been paid to bd transition-metal oxides, such as iridates 
(A2lr207), in search for unconventional phenomena that 
are induced by the competing spin-orbit (SO) coupling 
and electron correlation. 



Cd20s207 is one of the few compounds whose low-T 
magnetic structure has been experimentally determined 
among a number of magnetic bd pyrochlore oxides. The 
Cd^+ ion is non-magnetic, whereas the Os^+ ion with 
a bd^ configuration can be magnetic. This compound 
exhibits a purely electronic continuous metal-insulator 
transition (MIT) concurrently with a Neel ordering at 
Tmit ^ 227 K [i-lsl. Magnetic-susceptibility [1 % 
and /iSR measurements suggested the Neel ordering, 
whereas powder neutron diffraction did not confirm the 
abovementioned observation Q. Recently, Yamaura et 
at. have successfully detected a magnetic reflection at 
the wave vector of gr = using resonant X-ray scat- 
tering on high-quality single crystals f5|. They showed 
that only the so-called all-in/all-out magnetic order [see 
Fig. 1(a)] is compatible with the cubic symmetry of the 
crystal among the possible gr = magnetic orders. The 
all-in/all-out magnetic ordering has been suggested also 
in other bd pyrochlore oxides, i.e, Nd2lr207 in experi- 
mental work Q and Y2lr207 in theoretical work j9|. 

Despite the recent experimental progress, we are still 
far from fully understanding the electronic properties of 
this compound. In particular, several puzzling electronic 
properties have been reported experimentally: (1) In con- 
trast to the opening of an optical gap of the order of 
800 cm" -1 (ci 1100 K) [13, the semiconducting gap con- 
tinuously vanishes toward low T's [H, This indicates 
the absence of a clear charge gap at low T's which seem- 
ingly contradicts the semiconducting behavior of the re- 
sistivity up to Tmit- (2) The Neel transition temperature 
of this compound (= Tmit) is, to the beset of our knowl- 
edge, one of the highest among the magnetic pyrochlore 
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FIG. 1: (color online), (a) Cubic unit cell of Cd2 0s2 07 con- 
tains 8 formula units. Only the corner sharing tetrahedron 
network of Os atoms is shown for clarity. The arrows repre- 
sent the all-in/all-out magnetic order with the ordering vector 
of g = 0. (b)/(c) UeE dependences of the Os magnetic mo- 
ment mos, the charge gap Ac, and the direct gap Ad- The 
phase diagram consists of a non-magnetic metal (NMM), an 
antiferromagnetic metal (AFM) , and an antiferromagnetic in- 
sulator (AFI). 



oxides pLi]. This observation is quite surprising because 
geometrical frustration tends to prevent long-range mag- 
netic ordering. These puzzling characteristics and the 
similarity with the Ir oxides urged us to investigate this 
compound as a prototype of bd pyrochlore oxides. 

In this Letter, we investigate the ground-state proper- 
ties of Cd20s207 using extensive LSDA+S0+/7 calcula- 
tions (LSDA denotes the local spin density approxima- 
tion). We explore the ground-state phase diagram as a 
function of U^f^ [=17 — J). Here, [/eff is an empirical pa- 
rameter, and /7eff = 1-2 eV is expected to be appropriate 
for spatially extended 5d orbitals [9|]. We show that the 
all-in/all-out magnetic order is stable in an antiferromag- 
netic metal (AFM) and in an antiferromagnetic insula- 
tor (AFI). We observe that strong magnetic anisotropy 
arising from the SO coupling stabilizes the all-in/all-out 
magnetic order. Finally, we discuss possible origins of 
the experimentally-observed puzzling characteristics of 
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In the following calculations, we use a fully- 
relativistic two-component first-principles computational 
code, QMAS (Quantum MAterials Simulator) [pi. We 
employ the projector augmented-wave method [12[ and 
the LSDA-fSO+[/ method The relativistic effect 

including the SO coupling is fully considered in solving 
the relativistic Kohn-Sham equation 16, l^- In the fol- 
lowing calculations, we adopt a face-centered cubic primi- 
tive unit cell for Cd20s207 containing two formula units, 
and restrict consideration to q = magnetic ordering. 
Brillouin-zone integrations were performed using up to 
12x12x12 /c-point samplings using the improved tetrahe- 
dron method [l^ . We used a planewave cutoff energy of 
40 Ry. The following calculations were done with the ex- 
perimental lattice structure at 180 K: a = 10.1598 A and 
^(Oi) = 0.319 [3]. Every Os atom is located at the center 
of an OsOe octahedron. For x{Oi) > 0.3125, each oxygen 
octahedron is slightly compressed along the local (111) 
axis that connects the centers of the two neighboring Os 
tetrahedra. 

To identify the antiferromagnetic ordering, we calcu- 
late the local magnetic moment projected on an Os atom, 
mos5 by integrating the magnetic moment within a radius 
of 2.5 a.u. (= 1.323 A). The direct gap. Ad, is defined as 
the minimum gap between the conduction and valence 
bands, which approximately corresponds to the optical 
gap. In the following discussion, Ac denotes the charge 
gap, which identifies the MIT. 

Figure 1(b) shows the computed ground-state phase 
diagram with respect to /7eff. At small /7eff's, the ground 
state is non- magnetic metal (NMM). By increasing /7eff, 
the ground state turns into the AFM phase at /7eff 0.75 
eV, and further into the AFI phase at [/eff — 1-2 eV. 
We found that the all-in/all-out magnetic order is the 
most stable in the entire parameter region of the AFM 
and AFI phases: Os moments with the same magnitude 
point toward or away from the centers of the tetrahedra 
along the local (111) axes as illustrated in Fig. 1(a). We 
obtained mos — 0.8-1.1 jub/Os in the AFI phase, which 
is considerably smaller than 3 /ie/Os for the high-spin 
state. 

Figures 2(a)-(b) show the calculated electronic band 
structures near the Fermi level. At /7eff = 0.0 eV 
[Fig. 2(a)], the Fermi level lies in the half- filled t2g bands 
consisting of twelve Kramers (doubly )-degenerate bands. 
Note that the effect of the distortion of the octahedral 
crystal field is not strong enough to split the t2g manifold. 
As suggested by the previous LSDA studies [li,[2o|, the 
band structure looks semi-metallic on the high-symmetry 
lines. In fact, we obtained Ad — 0.035 eV at /7eff = 0.0 
eV. Although Ad decreases with increasing /7eff, Ad 
remains non zero within the NMM phase as shown in 
Fig. 1(c). Figures 3(a)-(f) show the calculated electronic 
density of states. At /7efr = 0.0 eV [Fig. 3(a)], one clearly 
see a sharp peak near the Fermi level, which may orig- 
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FIG. 2: (color online). Band structures on high-symmetry 
lines computed at UeS = 0.0, 0.5, 0.8, 0.9, 1.0, 1.1, 1.25, and 
1.5 eV. Energy, is measured from the Fermi level (^f) [(a)- 
(f)] or from the center of the charge gap [(g), (h)]. Conduction 
and valence bands are denoted by green and red lines, respec- 
tively. The direct gap closes near the L point at Z7eff — 0.9 
eV [see the circle in (d)]. 



inate in low- lying flat bands. The low-energy structure 
of the density of states remains essentially unchanged 
within the NMM phase [see Fig. 3(b)]. 

With increasing [/eff, the ground state turns into the 
AFM phase at Uef^ ^ 0.8 eV. We found that the all- 
in/all-out magnetic order is the most stable in the AFM 
and AFI phases. We used various initial spin and charge 
densities for the iterative scheme, but we could not find 
any other stable magnetic solutions. In the AFM phase, 
the time-reversal symmetry breaking lifts the Kramers 
degeneracy in the NM band structure. Note that the 
AF transition is associated with no Brillouin-zone fold- 
ing because the all-in/all-out order preserves the trans- 
lational symmetry of the lattice. The direct gap. Ad, 
closes at /7eff — 0.9 eV with a band inversion between 
the valence and conduction bands [see Fig. 1(c) and the 
circle in Fig. 2(d)] [33|. As mos develops with increasing 
/7eff, Ad opens again and becomes increasingly larger, 
finally resulting in the continuous MIT at [/gfr — 1.2 eV 
[Fig. 1(c)]. The MIT is characterized by vanishing elec- 
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FIG. 3: (color online). Density of states computed at [/gfr = 
0.0, 0.5, 0.9, 1.1, 1.25, and 1.5 eV for the same data as in 
Fig. 2. The inset of (e) is an enlarged plot of (e) near the 
charge gap. The density of states is measured per primitive 
unit cell. The shaded areas denote a pseudo gap, which ap- 
pears near the MIT. 



tron and hole Fermi surfaces (Lifshitz transition). This is 
clearly distinguished from the Slater transition in which 
a charge gap app ears because of a magnetic superlat- 
tice structure [21| . Another notable observation near the 
MIT is that the density of states is considerably sup- 
pressed near the Fermi level up to about 0.2 eV, which 
is remarkably higher than Ac [Figs. 3(d)-(f)]. As seen in 
Fig. 3(d), this pseudo gap starts to develop in the AFM 
phase, suggesting that it comes from the modification of 
the band structure near the Fermi level with the emer- 
gent magnetic order. 

Next, we discuss the magnetic properties of the AFI 
phase. It is well known that pyrochlore antiferromagnets 
tend to be frustrated when the spins are isotropic. This 
is clearly seen in the nearest-neighbor classical antiferro- 
magnet which exhibits no phase transition and remains 
paramagnetic down to zero T with a macroscopic degen- 
eracy [22|. The ground-state manifold consists of degen- 
erate states in each of which the summation of four spin 
moments vanishes on every tetrahedron. The degeneracy 
can be lifted by magnetic anisotropy which arises from 
the SO coupling. In particular, the all-in/all-out order 
is selected as the unique ground state by the local (111) 
easy-axis anisotropy [23h25i]. 

Figure 4(a) shows magnetic anisotropy energies cal- 
culated for Uef[ = 1.25 and 2.0 eV, respectively. Here 
Eg{0) is the energy of the self-consistent solution that is 



obtained under the constraint that every Os moment is 
rotated from that of the ground state around the [001] 
axis, and its magnitude is equal to mos in the ground 
state. One can clearly see that Eg remarkably increases 
with the rotation by as large as about 40 meV/Os and 
70 meV/Os for /7eff = 1.25 and 2.0 eV, respectively. This 
proves the existence of strong magnetic anisotropy in this 
compound even near the MIT. 

To gain deeper insight into the nature of the magnetic 
anisotropy, we extend the analysis using a phenomeno- 
logical model for the energy change associated with the 
rotations of the Os moments: 



U = J ^ rhi-mj - Asia ^ {rhi ■ dif 
+Adm ^ dij ■ {rfii xrfij), 



(1) 



where rhi is a unit direction vector of the magnetic mo- 
ment at Os site i. Here J, Asia, and Adm are the nearest- 
neighbor exchange interaction, the single-ion anisotropy, 
and the DM interaction, respectively. Note that mos 
weakly depends on /7eff in the AFI phase, and this effect 
is renormalized into the interaction parameters. The unit 
vector di is along the local (111) axis at site i. The case 
where Asia > corresponds to the easy-axis anisotropy. 
The unit vectors dij are the direction vectors of the DM 
interaction, which are the only one symmetry- allowed 
form [2^ Q . In Ref. the cases where Adm > and 
Adm < are distinctly referred to as the "direct case" 
and "indirect case"; the former favors the all-in/all-out 
ordering. 

The model parameters can be extracted by fitting en- 
ergies obtained in electronic structure calculations by 
Eq. ([1]). Hereafter, the energy is always measured per 
Os atom. In the following calculations, the magnitudes 
of the Os moments are fixed to mos of the ground state. 
On the basis of Eq. ([1]), the energy difference between 
the all-in/all-out and 3-in/l-out states is given by 



E^{{)) - Eg{{)) = J + 2a/2A 



DM- 



(2) 



Here, the 3-in/l-out state is the lowest excited state for 
J > and Asia > 0, which is obtained by flipping one of 
four non-equivalent spins in the all-in/all-out state [see 
Figs. 4(b) and (c)]. Changes in the energies of the all- 
in/all-out [EgiO)] and 3-in/l-out ordered states [Eq{0)] 
are obtained respectively as follows: 



^gW-^g(O) = ^eW-^e(O) 

4V^ 



(1 -cos6>). 



(3) 



EM -EM = Asia|l-^(2cos^ + l)'|, (4) 
where denotes the rotation angle around the [001] axis. 
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Table |I] summarizes the model parameters extracted 
for [/eff = 1.25 and 2.0 eV. The values of A^ia and Adm 
were obtained by fitting the magnetic anisotropy energies 
at smah (9's, i.e., in the range of 0° < 6> < 60° [35|. We 
found that the values of J and Asia are positive for both 
the abovementioned values of t/eff's. The antiferromag- 
netic J might originate in the superexchange coupling via 
the O site. The easy- axis anisotropy Asia is found to be 
even larger than J for both /7eff's. On the other hand, 
Abm is considerably smaller than Asia, and vanishes as 
C/gff increases. This might be because the DM interac- 
tion appears as a perturbation with respect to hopping 
between the nearest-neighboring Os atoms. These re- 
sults suggest that the all-in/all-out order is mainly stabi- 
lized by the large easy-axis anisotropy rather than the 
DM interaction. In the case of Y2lr207 [9], the all- 
in/all-out ordering is ascribed to the direct DM inter- 
action because Ir^+ has an effective total angular mo- 
mentum of Jeff = 1/2 [al, [23, [28|[36|. The observation of 
the large Asia indicates that Os^+ has a larger value of 

Jeff (Jeff > 1/2). 

Assuming that the real material is located near the 
MIT, e.g, f/eff ^ 1.25 eV, the calculated results provide 
a natural explanation for the peculiar low-T properties 
of this compound. (1) The large direct gap. Ad, and 
small charge gap, Ac, near the MIT are consistent with 
the experimental observation of a large optical gap and 
the absence of a clear charge gap at low T's. In fact. 
Ad — 0.17 eV estimated at [/eff = 1-25 eV is compara- 
ble to the optical gap of about 0.1 eV observed below 
Jmit 0- The pseudo gap is expected to develop be- 
low Tmit concurrently with the onset of the Neel order- 
ing. The density of states vanishes toward low T's in the 
wide energy range. This observation explains the semi- 
conducting behavior of resistivity up to the high Tmit- 
(2) The strong easy-axis anisotropy on the order of sev- 
eral tens meV stabilizes the all-in/all-out order cooper- 
atively with the antiferromagnetic exchange interaction. 
This can account for the high antiferromagnetic transi- 
tion temperature of this compound. The present result 
indicates that this compound has a magnetic gap of sev- 
eral tens meV. Further experiments are needed to detect 
the magnetic gap and confirm the present observation. 

Before closing this Letter, we would like to comment 
on future work. Contrary to the recent theoretical pro- 
posals 29I, 30| , all the three phases in the phase diagram 
have a trivial Z2 topological invariant. The band inver- 
sion in the AFM phase might suggest the presence of 
Dirac cones. A detailed analysis of the nature of the 
band inversion is a part of our future study. In con- 
trast to Cd20s207, Hg20s207 remains metallic below 
the AF transition temperature (sij . This material might 
be located in the AFM phase. The absence of a clear 
charge gap below the continuous MIT associated with the 
Neel ordering is also observed in NaOsOa perovskite [sij . 
Clearly, first-principle studies on this related compound 
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FIG. 4: (color online), (a) Magnetic anisotropy energies es- 
timated by the rotation of Os magnetic moments around the 
[001] axis. (b)/(c) all-in/all-out and 3-in/l-out ordered states 
rotated by 6* = 0° and = 120° around the [001] axis. Every 
spin is perpendicular to the local (111) axis at ^ = 120°. In 
(b), the small (gray) vectors and the labels represent the di- 
rection vector 
i, respectively. 
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TABLE I: Estimated parameters for the phenomenological 
model given in Eq. ([T]). 



will be interesting. 

In summary, we performed LSDA+SO+[/ calculations 
to explore the ground-state properties of Cd20s207. We 
also found the all-in/all-out magnetic order as the ground 
state in a wide region of the phase diagram, support- 
ing recent X-ray experimental results. We also showed 
that this magnetic order is stabilized by strong easy-axis 
anisotropy originating in spin-orbit coupling. Further- 
more, we found a pseudo gap in the density of states 
near the metal-insulator transition. These numerical re- 
sults provide a natural explanation for the puzzling low- 
temperature properties of this compound. The present 
result might open up new possibilities to explore Os com- 
pounds as a stimulating playground for the interplay of 
electron correlation and spin-orbit coupling. 
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